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ABSTRACT 

Context. Some O stars are suspected to have to have (weak) magnetic fields because of the observed cyclical variability in their UV 
wind-lines. However, direct detections of these magnetic fields using optical spectropolarimetry have proven to be very difficult. 
Aims. Non-thermal radio emission in these objects would most likely be due to synchrotron radiation. As a magnetic field is required 
for the production of synchrotron radiation, this would be strong evidence for the presence of a magnetic field. Such non-thermal 
emission has already been observed from the strongly magnetic Ap/Bp stars. 

Methods. We have performed 6 & 21 cm observations using the WSRT and use these, in combination with archival VLA data at 3.6 
cm and results from the literature, to study the radio emission of 5 selected candidate magnetic O stars. 

Results. Out of our five targets, we have detected three: £ Per, which shows a non-thermal radio spectrum, and a Cam and A Cep, 
which show no evidence of a non-thermal spectrum. In general we find that the observed free-free (thermal) flux of the stellar wind 
is lower than expected. This is in agreement with recent findings that the mass-loss rates from O stars as derived from the Ho- line are 
overestimated because of clumping in the inner part of the stellar wind. 
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All O-type stars have strong, line-driven winds. They usually 
have a thermal radio spe ctrum due to the free-free emission from 
the ionised stellar w ind (lAbbott et al.lll980tlBieging et al.ll 19891: 
IScuderi et alJll998l ). This spectrum can be calculated using: 
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(IWright & Barlowl [1975 iPanagia & FeiH [1975 IScuderi et alJ 
1998j), where D is the distance to the star, M is the mass-loss 
rate, T e the electron temperature, yU e the mean atomic weight per 
electron, v«, the terminal wind velocity and v the observing fre- 
quency. 

However, about 30% of the O s tars are found to show non- 
thermal radio emission (see, e.g.. |Bieging et al] 119891 : iDrakd 
Il990t IScuderi et alJll998t iBenaglia et al.ll200lh . This is charac- 
terised by a flat ter than therm al spectrum, i.e. defined as a < 0.6, 
with S v oc V. I White! ([T985) proposed that synchrotron radia- 
tion from the rapidly moving electrons of the wind in a stel- 
lar magnetic field could also contribute to the radio emission. 
This was confirmed by the discover y of non-thermal r adio emis- 
sion in the mag netic Ap/Bp stars (iDrake et alJll987l see also 
ICassineliilll984l who report non-thermal emission found in <x 
Ori E by Churchwell). These stars have strong (of the order of 
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kilogauss), dipole-like magnetic fiel ds. Among O stars only two 
magneti c stars are known: 6 l Ori C (iDonati et al.ll2002b and HD 
191612 dDonatiet all 12006). Nevertheless, strong indirect evi- 
dence exists that ma ny O stars should have magnetic fields (e.g. 
iHenrichs et al.ll2005) . One of the main arguments why many O 
stars are thought to have (weak) magnetic fields is that their 
winds show cyclic beha viour on a rotat ional timescale, which is 
typically a few days (see lFullertonl2003l for a review). The lack 
of magnetic field detections is most likely related to the fact that 
direct measurements of magnetic fields in O stars are extremely 
difficult, because of the very few available spectral lines in the 
optical region. The usual method to measure magnetic fields is to 
determine the magnetic Zeeman splitting of magnetically sensi- 
tive lines with optical spectropolarimetry. The sensitivity of this 
method decreases towards earlier spectral types, as these stars 
have fewer spectral lines in the optical region. 

The detection of non-thermal radio emission from O stars 
with such cyclic variability would be strong evidence that mag- 
netic fields are indeed present in these stars. We have selected 
five candidate magnetic O stars that have been studied exten- 
sively in the ultraviolet (UV) in order to search for evidence of 
non-thermal radio emission. 



2. Observations & data reduction 

For this study five targets have been selected from the 10 O 
stars listed by Kan er et all (Il996l) . which are the brightest and 
best studied O stars in the UV with the International Ultraviolet 
Explorer (IUE) satellite. All these targets show extensive stel- 
lar wind variability, some with well studied cyclic behaviour. 
The final selection was made on the criteria that the star should 
be observable with the Westerbork Synthesis Radio Telescope 
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Tabl e 1. Stellar parameters of the selected O stars. Spectral types are from Walborr] dl973Lfl976l) . Hipparcos parallaxes were taken 
from Perrvma n et al.l (Il997l) . All other paramet ers were taken from th e "preferred solution" (some stars have more than one possibl e 
set of parameters that fits the observations) of iMarkova et al.l (120041) . except for 10 Lac, for which we used lMokiem et alj (2005). 
We show both the distance from Hipparcos and the distance used for the spectral modeling. As this latter value is us ed to determine 
the ma ss-loss rates, this distance was also used for predicting the thermal radio flux. For the preferred solutions of Mar kova et al.l 
(12001 . we have scaled the distance from the original solution using the absolute magnitudes, as the reddening is assumed to be the 
same for both solutions. 





£Per 


a Cam 


15 Mon 


/ICep 


10 Lac 


HD number 


24912 


30614 


47839 


210839 


214680 


Association/Runaway 


Runaway 


Runaway 


Mon OBI 


Runaway 


Lac OBI 


Spectral type 


07.5III(n)((f)) a 


09.5 la 


07 V ((f)) 


06 I (n)fp 


09 V 


Parallax (mas) 


1.84 ± 0.70 


0.47 ± 0.60 


3.09 ± 0.53 


1.98 ±0.46 


3.08 ± 0.62 


Hip. distance (pc) 


540+ 330 


>440 


323 + ^ 


510 +15 ° 


320+* 


Spectral mod. distance (pc) 


895 


765 


710 


1089 


320 


Mass (Mq) 


52 


22 


25 


58 


27 


Radius (R Q ) 


25.2 


19.6 


9.9 


23 


8.3 


7-eff(10 3 K) 


34.0 


31.0 


37.5 


36.2 


36.0 


Luminosity (10 5 L o ) 


7.6 


3.2 


1.7 


8.1 


1.0 


Terminal velocity (km s ) 


2400 


1550 


2200 


2200 


1140 


Mass-loss (1(T 6 M y" 1 ) 


4.0 ± 1.0 


2.9 ± 0.9 


1.2 ±0.3 


7.7 ±2.3 


6.1+^ x 10- 2 



This is the spectral type given by Walborn ( 1973), however, a spectral type of 07.5I(n)((f)) was adopted by Mar kova et all {2004). 



(WSRT) and that the star should have been previously detected 
in the radio region (a Cam, 15 Mon and A Cep; for £ Per archive 
observations from the Very Large Array -VLA- were available). 
We added 10 Lac because of its brightness and its rich UV his- 
tory. The stellar parameters of these stars are listed in TableQ] 

The radio observations of our targets selected from the lit- 
erature have been summarised in Table [2] The detection of 
a 42 + 5 mJy source nea r the o ptical position of £ Per by 
iBohnenstengel & Wendkerl d 19761) at 1 1 cm (2.7 GHz) is dis- 
cussed in Sect. 13. II 

To complement these measurements, and in order to deter- 
mine the spectral slopes, we have used WSRT (6 and 21 cm) and 
VLA (3.6 cm) observations. 



Table 2. Radio detections as a function of wavelength of our 
targets reported in the literature. When several measurements 
are available a weighte d average is shown ; up per limits are 
2,cr . Fl uxe s are from [11 lAbbott et all d!980l). [21 [Bieging etal] 
(I1989I). [3llDrake1(ll990l) . [4l lLamers & Leithererl(ll993l) and 151 
IScuderietalJ dl998). 



Star 


Flux (mJy) 






2 cm 3.6 cm 6 cm 


References 


a Cam 


0.65 ±0.13 0.44 ±0.04 0.29 ± 0.04 


2,5 


15 Mon 


<0.4 0.40 ±0.13 


3 




<0.33 & <0.18 


2.3 


A Cep 


0.38 ± 0.03 0.40 ± 0.25 


1,4 



We performed observations for all five selected O stars at 21 
cm (1.4 GHz) with the WSRT during the period from September 
to November 2005. In addition, 10 Lac was observed at 6 cm 
(4.9 GHz). All observations consisted of 12 h integrations in the 
Maxi-Short configuration, done in continuum mode with a band- 
width of 8x20 MHz. Gain and phase calibrations were done us- 
ing the calibrator 3C286, except the observation of a Cam which 
was calibrated with 3C48. 



Earlier observations of £ Per were performed in 1995 at 6 cm 
(6 observations in May and June, total of 26 h) , and 21 cm (5 
observations in June, July and August, total of 39 hours). Due to 
the lower sensitivity of the WSRT at that time, the observations 
at each frequency were all combined. The 21 cm observations 
were calibrated using 3C48 and the 6 cm observations were cal- 
ibrated using 3C48, 3C147 and 3C286. 

The reduction of the WSRT data was done using the miriad 
software package. 

From the VLA archive we used an X-band (3.6 cm, 8.5 
GHz) continuum observation taken on 11 Jan 1999 (program 
ID AS644-x) of £ Per. This 0.68h observation was taken in C 
configuration with a bandwidth of 50 MHz. The data were re- 
duced with aips, using 3C48 as a primary and B041 1+341 as a 
secondary calibrator. 

2.1. Distances and mass-loss rates 

As massive stars are relatively far away, their distances as de- 
termined by the Hipparcos satellite suff er from systematic er - 
rors and and may be underestimated (e.g. Schr oder et al.l l2004). 
Distances may also be obtained using the stellar properties de- 
rived from spectral modeling. The spectral modeling distance for 
a star is generally derived from the relation between the spec- 
tral type and luminosity, possibly making use of other stars in 
the same cluster. This distance is used to determine the absolute 
magnitude of the star and this constrains the stellar radius. As the 
mass-loss rates determined in this way depend on the radius, we 
have adopted the distances used in the spectral modeling to cal- 
culate the predicted radio fluxes. This ensures that our mass-loss 
rates and distances are consistent (see also Sect.©. 



3. Results 

We detected three of our five selected targets (see Table [3}. For 
£ Per this was the first detection in the radio; all three stars have 
now been detected for the first time at 2 1 cm. In general the flux 
was found to be lower than the predicted thermal flux using Eq.Q] 
and the stellar parameters from Table[T](assuming typical values 
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of He = 1.3 and T e = 0.85 r, ff : IScuderi et~ai1ll998l) . Both the 
predicted thermal flux and a power law fit to the observations 
are shown in Fig.Q] 

We now present the results per source in order of RA: 

3.1. £ Per 

In the 1995 WSRT observations, £ Per was not detected at 6 
and 21 cm. However, it was detected in the higher S/N 3.6 cm 
(VLA, Jan 1999) and 21 cm (WSRT, Nov 2005) observations. 
The flux was found to be lower than the predicted thermal flux 
by a factor of ~2 (21 cm) to ~3.5 (3.6 cm). The spectrum has a 
spectral index of a — 0.29+0. 14, which is lower than the thermal 
value of 0.6. This is evidence for the presence of a non-thermal 
contribution to the ob served flux. 

iPuls et alJ d2006l) found an upper limit for £ Per of 120 //Jy 
(3cf) from VLA observations on March 9, 2004. As this limit is 
not consistent with our detection at this wavelength, this might 
be an indication of variability. To check this, we retrieved the 
observations from the VLA archive. At the position of f Per, we 
measured a flux density of 154+39 //Jy, which, we agree, is not a 
reliable detection of the source (3.9<x). However, it is consistent 
with o ur detection of the source at 1 69 ± 30 //Jy in 1999. 

IBohnenstengel & Wendkerl (fl976) detected a source of 42+5 
mJy at 2.695 GHz near the optical position of £ Per. They con- 
cluded that this component is either due to an extended (~2') 
thermal source of about 10 mJy, or to blending of their compo- 
nents A and B (see Fig. 0, in which case they claim that com- 
ponent B has to have a very flat spectrum. As an interferometer 
such as the WSRT is not very sensitive to extended structures, 
we cannot exclude the presence of an extended source. We find 
that component B has a spectral index of orw-1.1. It is detected 
at 5.8+0.2 mJy at 21 cm, its 6 cm flux is 1.5+0.1 mJy and at 
3.6 cm it shows extended structure (10"x2") but has very low 
flux density (peak of ~0.2 mJy/beam). At 3.6 and 6 cm, compo- 
nent A can be resolved into two components with a separation of 
12+1". The western component hasafluxof6.1±l.l(6cm) and 
2.6+0.5 mJy (3.6 cm) and the eastern component 7.1 ±0.9 (6 cm) 
and 3.0 ± 0.5 mJy (3.6 cm). In addition we found a third source 
at the position [a(2000)=03h59m03s, 5(2000)=+35°46'42"], 
which is not detected at 6 cm (< 0.24 mJy) but is detected at 
0.9+0.2 mJy at 21 cm. 

Given the low resolution of the Effelsberg telescope (~5' at 
1 1 cm) compared to the WSRT (-15" at 21 cm) and VLA (~2" 
at 3.6 cm) observat ions, we conclude that it is very l ikely that 
the source found by iBohnenstengel & Wendkerl (1 19761) is due to 
blending of all these components. 



xi Per 21 cm 




3 h 59 m 20= 59 m 10 s 59 m 00 s 58 m 50= 5B m 40 ! 



Right Ascension (J2000) 

Fig. 2. Our 21 cm $ Per observation of 2005. The sources A and 
B of iBohnenst engel & Wendkerl d 19761) are indicated, as are £ 
Per and a new source that they did not detect (see Section [3~TT >. 



Due to the unfavourable position, for an E-W array, of 15 
Mon on the sky (close to the equator), our 5<x-upper limit of 250 
//Jy at 21 cm is relatively high compared to the expected thermal 
noise for a 12h run. 

3.4. A Cep 

Our detection of A Cep at 21 cm allows for an accurate de- 
termination of the spectral index. We find that A Cep has an 
approximately thermal spectrum, with a spectral index of a = 
0.74 ± 0.1 1. The flux is, however, a factor of -3-3.5 lower than 
predicted. 

3.5. 10 Lac 

We have not detected 10 Lac, with a 5cr upper limits of 95 
//Jy at both 21 and 6 cm. This is in agreement with the ex- 
pected ther mal flux based on th e determination of the mass- 
loss rate bv lMokiem et all d2005l) of 6.1^| | x IO^MqV -1 . Our 
6 cm upper limit constrains the mass-loss rate to be lower than 
1.3 x 10- 7 M Q y _1 . 



3.2. a Cam 

This star shows a thermal spectrum over the entire range from 2 
to 21 cm. The flux is found to be lower than the predicted thermal 
flux by a factor of ~2. 

3.3. 15 Mon 

IDrakd (1 1990I) found a 3<x upper li mit of 400 //Jy on 24 Jan 1987 
for this source at 2 cm. At 6 cm Bieg ing et al.l dl989l) reported 
a 3cr upper limit of 330 /Jy, and iDraka (1 1 9901) found an upper 
limit of 180/Jy on 22 Feb 1986. However, one year later (24 Jan 
1987)|Drake| d 1990b detected the source at the same wavelength 
at 400 + 130 //Jy. This suggests that the radio flux of 15 Mon is 
possibly variable. 



4. The effects of clumping in stellar winds 

As re cently discussed by iFullerton et alJ (|2006) and IPuls et al.l 
(2006), it is generally found that mass-loss rates determined 
from Ha are higher by a factor of ~3-8 than those derived from 
radio observations. This is thought to be due to enhanced clump- 
ing in the inner part of the wind where the Ha emission origi- 
nates, compared to the outer part of the wind from which we 
receive the observed radio emission. As both the Ha and radio 
emission are proportional to the density squared, clumping in the 
wind results in an overestimate of the mass-loss rate. Since we 
use mass-loss rates derived from Ha to predict the thermal radio 
fluxes, the fact that the observed fluxes are lower than our pre- 
dictions is consistent with stronger clumping in the inner part of 
the wind compared to the outer part. 
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Table 3. Results of our new WSRT and archival VLA observations. The upper limits shown are 5cr upper limits. 



HD 


star 


date 


freq. 


A 


flux 


spectral 


array 


number 


name 


(d/m/y) 


(GHz) 


(cm) 


fiJy 


index a 




24912 


fPer 


May/Jun/95 


4.9 


6 


<240 




WSRT 






Jun/Aug/95 


1.4 


21 


<190 




WSRT 






ll/Jan/99 


8.4 


3.6 


169+30 


0.29+0.14 


VLA 






28/Nov/05 


1.4 


21 


100+18 




WSRT 


30614 


a Cam 


09/Oct/05 


1.4 


21 


156±12 


0.57±0.06 


WSRT 


47839 


15 Mon 


13/Oct/05 


1.4 


21 


<250 




WSRT 


210839 


A Cep 


03/Oct/05 


1.4 


21 


98±21 


0.74±0.11 


WSRT 


214680 


10 Lac 


22/Sep/05 


4.9 


6 


<95 




WSRT 






21/Sep/05 


1.4 


21 


<95 




WSRT 



!-3 




15 Mon 




07 V ((f)) 


V 


V 









0.60.8 1 2 4 6 8 10 20 — 0-60.8 1 2 4 6 8 10 20 0.60.8 1 



3 






E i 


X Cep 






61 (n)fp 




0.1 







4 6 8 10 20 




• new detection 

* literature 

v upper limit (5a) 



4 6 8 10 20 



4 6 8 10 20 



Frequency (GHz) 

Fig. 1. Radio spectra of the 5 selected targets. Shown are the new results, results taken f rom the literatur e, upper limits, the predicted 
thermal flux (Eq. [TJ dashed line), the predicted thermal flux accounting for clumping dPuls et alJl2006l dotted line), and a fit to the 
observations (solid line). Distances used for estimating the thermal flux are the same as those used for the determination of the 
mass-loss rates using Ha. 



We have used the same distances as assumed in the Ha mod- 
eling. For the Ha modeling the distance (or absolute magnitude) 
is used to estimate the stellar radius, resulting in R oc D. In 
these models the mass-loss rate approximately scales with R as 
M oc R?l 2 , which gives M oc D 3 ^ 2 . Since the predicted radio 
flux scales as S v oc M 4 / 3 D~ 2 (Eq.[TJ, one finds that the predicted 
radio flux is approximately independent of the distance. 

For a Cam and A Cep, our results ar e in goo d agre ement 
with the mass-loss rates determined by iPuls et all (|2006), and 
for £ Per, the flux observed at 3.6 cm is in agreement with the 
upper limits of both possible solutions quoted. 

5. Conclusions 

We have detected three candidate magnetic O stars at radio 
wavelengths. Of these three £ Per shows a non-thermal spectrum 
and a Cam and A Cep show thermal spectra. As non-thermal ra- 
dio emission is assumed to be due to synchrotron emission, the 



detection of a non-thermal radio spectrum in £ Per strengthens 
the case that the observed UV line variability observed in this 
star is caused by a magnetic field. 

5.1. The origin of non-thermal radio emission 

IVan Loo etail d2004 see also iRauw et al] I2002I) were able to 
successfully reproduce the non-thermal emission of the 05 If 
star Cyg OB2 No. 9, using a model assuming a power-law for the 
momentum distribution of the relativistic electrons. Reasonable 
values for the magnetic field that were used are of the order of 
10-100 gauss. However, more recent numerical simulations in- 
corp orating individual sh o cks for the acc eleration of electrons 
(e.g. Ivan Loo et a n i2005t IvanLool 12005) suggest that both a 
magnetic field and a binary companion are required to explain 
non-thermal radio emission from massive stars. In these simu- 
lations, the synchrotron radiation from single massive stars with 
magnetic fields is produced relatively close to the star, where 
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shocks occur which accelerate the electrons. This radiation is 
absorbed in the stellar wind due to the large free-free opac- 
ity. When a star has a massive binary companion, electrons are 
accelerated in the wind-wind collision region, where the radio 
emission can escape because of the lower opacity (r, a( jio < 1 ). 

In the less massive strongly magnetic Bp stars, non-thermal 
radio emission is observed that is thought to be produced by 
mildly relativisti c electrons which are trapped in the stel lar mag- 
netosphere (e.g. brake et al.1 Il987t iLeone et al.1 Il996l) . Radio 
observations of known magnetic massive stars would help to 
increase our understanding of the relation between magnetic 
field strength, mass-loss rate and non-thermal radio emission. 
Unfortunately, the few massive stars that have measured mag- 
netic fields are likely relatively weak radio sources as they have 
rather weak (order 10 2 gauss) fields and a low M (J3 Cep, V2052 
Oph, ( Cas and r Sco), or are in a crowded field (8 l Ori C). 

5.2. £ Per, A Cep and a Cam 

£ Per is a single runaway star dGies & Boltonl lT986). Runaway 
stars are stars with high space velocities of ~30-200kms~ 1 and 
are usually single stars. Therefore, at least for some single mas- 
sive stars it is possible to have a (mildly) non-thermal radio spec- 
trum. The detection of variability in the radio flux of close binary 
stars suggest th at stellar winds are not as optically thick as gen- 
erally assumed (Blomme 2005). Due to porosity effects, clump- 
ing and asphericity of the mass loss due to magnetic fields, it 
might be possible to observe radio emission from much closer 
to the central star. For producing synchrotron emission, a mag- 
netic field is necessary but not sufficient, as relativistic electrons 
are also required. In the case of £ Per, and other possibly single 
massive stars with non-thermal radio spectra, the precise origin 
of the relativistic electrons needs further investigation. One pos- 
sible origin in the case of runaway stars would be the bow shock 
of the star moving through the interstellar material. However, in 
high resolution radio observations, such as our 3.6 cm observa- 
tion of £ Per, such a bow shock would most likely be resolved. 

In principle, variability of the radio flux of £ Per could 
change the spectral index since the data at different wavelengths 
are taken at different epochs. However, observed radio variabil- 
ity of massive stars is usually related to the orbit of a massive 
companion. Such variability is not expected for £ Per, but as the 
mechanism responsible for the relativistic electrons is not com- 
pletely understood we plan future observations to confirm the 
non-thermal character of the radio spectrum, and check if vari- 
ability is present. 

Th e runaway stars A Cep and a Cam, which are presumably 
single dGies & Boltonlfl986i) . have a thermal radio spectrum, but 
it can not be excluded that they have a magnetic field. A Cep has 
a much denser stellar wind, and the non-thermal emission might 
all be absorbed. In the case of a Cam, the magnetic field might 
be too weak to produce observable non-thermal emission, but the 
lack of non-thermal emission could also be due to the location 
(closer to the star) or strength of the shocks required to produce 
the relativistic electrons. 

Finally, o ur results agree with r ecent results by 

iFullerton et al.l d2006h and iPuls et all d2006h that the mass-loss 
rates as derived from free-free radio emission are significantly 
lower than those derived from Ha modeling, which is a 
signature of enhanced clumping in the inner part of stellar wind. 
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